The axial breathing mode in rapidly rotating Bose-Einstein condensates and 
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Experiments on the axial breathing mode in a rapidly rotating Bose-Einstein condensate are ex- 
amined. Assuming a cold cloud without thermal component, we show that errors due to defocus of 
an imaging camera in addition to an inclination of the rotational axis can lead to a significant un- 
derestimate of the rotation rate in the fast rotation limit; within these uncertainties, our theoretical 
prediction agrees with the experimental data. We also show that, in the fast rotation regime, the 
Thomas-Fermi theory, which is inapplicable there, underestimates the rotation rate. Underestimate 
of the rotation rate due to these effects would also partly explain a discrepancy between theory and 
experiment for the Tkachenko mode frequency in the fast rotation regime. 
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PACS numbers: 03.75.Kk, 05.30.Jp, 67.40.Vs, 47.37.+q 



The creation of vortices in atomic Bose-Einstein con- 
densates has opened up a new horizon in the study of 
superfluids. This system enables us to investigate an un- 
explored regime, where the vortex core is comparable to 
the intervortex distance since the interaction energy per 
particle ~ gn is much smaller than that of liquid He-II 
(see, e.g., Ref. [l[ and references therein). Here n is the 
particle number density, g = 4iTh 2 a s /m is the two-body 
interaction strength, m is the particle mass, and a s is the 
s-wave scattering length. 

For a harmonically trapped Bose-Einstein condensate, 
when the rotation angular velocity f2 is close to the trans- 
verse trap frequency u>±, so that TiVL > gn, the condensate 
wave function is dominated by the lowest Landau level 
(LLL) component @. When the number of particles N 
is much larger than the number of vortices N v , the sys- 
tem may be described by the Gross-Pitaevskii equation 
[3, H[. Recently, Schweikhard et al. [f| have reached 
what we shall refer to as the mean-field LLL regime 
0, @, 0, B S E3, [LI G3, d, in which Hfl > gn with 
N>N V . 

In Ref. [f|, Schweikhard et al. have measured the ax- 
ial breathing mode frequency and they have observed a 
frequency shift in the rapid rotation regime, which can- 
not be described by the superfluid hydrodynamic models 
of Refs. [3, We have studied breathing modes in 
the mean- field LLL regime in Ref. [l6[ , and have pointed 
out that our theoretical framework developed there could 
describe this frequency shift. The purpose of this paper 
is to compare in detail our theoretical prediction of the 
axial breathing mode frequency with experimental data. 

In the Appendix of Ref. [10], in which we discussed 
the breathing modes in three dimensions, we employed a 
generalized LLL wave function ^ ex , in which the trans- 
verse oscillator length is treated as a variable. There 
we assumed the coarse-grained density profile n(r) = 
N(\4> ex \ 2 } to be a Thomas-Fermi (TF) parabola in the 
transverse direction and a Gaussian form in the axial (z) 
direction, which is correct in the fast rotation limit, in 
which the cloud is so dilute that huj z 3> gn, where lj z is 



the axial trap frequency. In the present work, we employ 
an extended version of this framework, which describes 
the coarse-grained density profile of the cloud more cor- 
rectly at lower rotation velocities. Here we assume 



n(r) = n(0) 1 



1R 2 



(1) 



where r 2 , — 



y 2 , R± is the TF radius in the trans- 



verse direction, a is a parameter (static), which deter- 
mines the shape of the axial density profile, and \faR z 
is the cloud radius in the axial direction. To satisfy 
the normalization condition, J d 3 r n(r) — N, we set 
ra(0) = 2N/[%y/aB (a + 1, ±) R 2 ± R Z ], where B is the 
beta-function defined as B(p,q) = f dt t p_1 (l — t)^ 1 . 
The axial density profile of Eq. (Q]) interpolates between 
the one-dimensional TF parabola (a = 1) and the Gaus- 
sian form (a = 00). 

The extended LLL wave function correspond- 
ing to the density profile |T]) is </> C x(r) = 
^e x nt v i(C - CO exp [- (A- 2 - %0) ri/24] [1 - 

{z 2 /aR 2 )\ a ^ 2 exp [i^z 2 /2d 2 z ] , where ( = x + iy, Q are 
the vortex positions, A cx is the normalization constant, 
d± = s/h/mux and d z = y/fi/muj z are the transverse 
and the axial oscillator lengths, respectively. Here the 
dynamical variable A describes the variation of the 
effective transverse oscillator length. The dynamical 
variables /3 and 7 generate a radial and an axial velocity 
field, respectively, which cause the homologous change 
of the density profile. 

Following the same way as in Ref. [l6| , we obtain the 
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Lagrangian function for the above <^ ex as 




+ ( 7 2 + l) h(a) 



where Xq and Zq are the equilibrium values of X and Z, 
and f2o is the rotation velocity in the equilibrium state, 
which we obtain as flo/uj± = 31 /Xq. Eliminating k from 
Eqs. (0) and (01, we write Qq as 



fi 



J 1 - 3 — w(M«)^o 4 -g(«)) • (5) 



(the dot denotes the time derivative) with /(a) = 
v / 2 7 r/a J B(2a + 1,±)/B 2 (a + 1,|), g(a) = 2aB(a - 
l,\)/B{a + /i(a) = 2aB(a + 1,§)/B(a + 

1, i), and the dimensionless interaction strength k = 
(v / 2^d z )- 1 (m6 ff iV)/(2^ 2 ), where 6 = (|<M 4 )/<l</>c X | 2 ) 2 
is the Abrikosov parameter, which is comparable to unity. 
Here X = R±/d±, Z = R z /d z , and I is defined by the 
expectation value l z of the angular momentum in the z 
direction per particle as / = {l z /K) + 1 — X 2 /3A 2 . Note 
that, in the limit of the axial density profile being Gaus- 
sian (a — > oo), the functions /(a), g(a), and h(a) are 
unity and thus Eq. © reduces to Eq. (A8) in Ref. [16| . 

Using the variational Lagrangian formalism, we obtain 
coupled equations similar to Eqs. (A. 12) and (A. 13) in 
Ref. [lj| that determine R± and R z (or a) in the equi- 
librium state: 



On 



1 - -r ) x o ~ 8«/(a) V = o , 



h(a)Z 3 



(a)Z n 



8k ft \ u± y- 

O U) z 



(3) 
(4) 



In addition to Eqs. ^ and Igj), we get an extra equa- 
tion by minimizing the energy with respect to a: 



h'{a)Zl+^{a)Z^ + ^ /'(«)— A ~ 2 = , (6) 



where the prime denotes the derivative with respect to 
a. Using Eqs. (T5|) and ©, Z is written as a function of 
only a: 



, , 2 <S(<*) lg!(a)y /A ( 2 h(a) , lh'(a)Y 1/4 



/(«) 2 /'(a) 



/(«) 2 /'(a) 



(7) 



Resulting expression of the mode frequency is 



± 



n 2 



4w 2 



+ 8w 2 u; 2 1 



g(a) 1 ' 
M«)^o 4 . 
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The 1/ Zq terms come from the zero-point energy in the z 
direction, which is not taken into account in the hydrody- 
namic calculations [3, UH based on the TF theory [ijj . 
We see that, in the limit of a — ► oo, Eq. ([8]) reproduces 
Eq. (A16) in Ref. [l6j], which is obtained for a cloud with 
a Gaussian form of the axial density profile [l8[. Note 
that Eq. ([SJ) can describe the frequency shift of the axial 
breathing mode from ySui z to 2ui z with increasing f^o in 
the fast rotation regime reported in Ref. [5[ . For values of 
f^o at which the frequency shift occurs, the second term 
in the square root of Eq. © can be neglected and this 
equation reduces to 

Wt 2 r ^i and o4=i(3 +f||^r) ^ W 

where u>t T and w ax are the transverse and axial breath- 
ing mode frequency, respectively. Since Z — R z /d z > 1, 



the first case of Eq. © is given by the upper sign of 
Eq. (JSJ) and the second one by the negative sign for 
uj± > lu z as in the experiment of Ref. § . If the system 
is in the LLL limit, z dependence of the wave function 
corresponds to the ground state of a particle in a har- 
monic potential, which is reproduced by a — > oo giving 
n(r) oc exp (— z 2 /i? 2 ) with R z = d z , then Zq ~ 1, and 
g(oo) = h(oo) = 1. Thus the second expression of Eq. 
© leads to w ax ~ 2uj z . Otherwise 1/Zq term is small 
and this expression gives w ax ~ V3ui z . This observation 
shows that the zero-point energy in the z direction, which 
is neglected in the TF theory, plays an essential role in 
the frequency shift of the axial breathing mode. We also 
note that, from the first case of Eq. ©, such a change 
in the mode fre quen cy does not exist for the transverse 
breathing mode [19| . 

In comparing experimental data and theoretical pre- 
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diction, a key issue is the determination of the rotation 
velocity f2 - In Ref. 0] , S7 is evaluated from the aspect 
ratio of the cloud using the TF theory, which is invalid in 
the fast rotation regime. Thus we do not use the data of 
Slo given in Ref. [5j; instead, we determine f2o using the 
present theoretical framework. In the experiment of Ref. 
[5j, what Schweikhard et al. have observed is the phase 
contrast image of the cloud, which is proportional to the 
optical density n op t(y,z) = J dx n(r) [20(. To perform 
direct comparison with the present theory, Schweikhard 
reanalyzed the stored data of the image obtained in that 
experiment and newly measured the rms radii of the op- 
tical density profile [201 ] . 

The rms radii are extracted in the following way. Since 
the signal to noise ratio is low in the fast rotation regime 
(see, e.g., Fig. 1(c) in Ref. 5]), contribution from the 
outer region of the cloud will produce a large error of the 
rms radii. To make the measured values immune against 
the noise in the images, image pixels with an optical den- 
sity lower than 50 % of the peak optical density had to 
be discarded [2(J. Thus the rms radii (y 2 ) 1 ^ 2 and (z 2 ) 1 / 2 
are measured for a part of the cloud in which the optical 
density is higher than 50 % of its peak value [20], i.e., 
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J K dydz y 2 n opt {y,z) 
J K dydz n opt (y,z) 



(10) 



where the integration region 1Z = {{y,z) | n opt (y,z) > 
0.5n op t(0, 0)} ((z 2 ) is given by an analogous expression). 

On the other hand, a theoretical expression of the op- 
tical density for the density profile of Eq. |T]) is 



(■opt 



(y,z) = -n(0)R x 1- 
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FIG. 1: (Color online) Axial breathing mode frequency cj ax 
as a function of fio- The horizontal error bars are the ex- 
perimental data and the points show our theoretical results 
(see text for details). The dashed line shows the result of 
hydrodynamic models by Refs. [3, [IH- 



the imaging resolution gets worse by a factor of y2 at 
most [2fJ. When the focusing is exact, an effect of the 
finite imaging resolution is corrected by subtracting in 
quadrature the rms width of the Gaussian point spread 
function (PSF) of 2.4 /mi from the raw data of the rms 
radii. Defocusing by one depth of focus corresponds to a 
point spread function of 3.4 /jm instead of 2.4 fim. [20| . 
Since (z 2 ) ex <C (y 2 )cx at f^o on which we are focusing, 
the ratio of increase of {z 2 )]J 2 due to the defocusing is 
larger than that of (y 2 )li 2 . Thus this effect makes an 
observed image more prolate than reality and leads to an 
underestimate of Qq. 

The other source is an overestimate of ( Z )ex ' due to a 
tilt of the rotation axis towards/away from the imaging 
camera. The tilt angle in the direction of the camera can- 
not be measured and thus the uncertainty of the tilt angle 



(11) 

Similar to the above experimental procedure, theoretical 
values of the rms radii are calculated from Eq. (|10p with 
Eq. ([ITp. We note that since the equilibrium values of 

R z and a are related by Eq. ©, (y 2 ) andjz 2 ) given by produces an error of (^1/2 _ Judging from the observ ed 

inclination in the plane perpendicular to the camera, the 
tilt towards/away from the camera is estimated at up to 
~ 5° [2^. This effect also leads to an underestimate of 
f2o due to the same reason as mentioned above. 

In Fig. [TJ we plot the resulting axial breathing mode 
frequency function of Qq. The horizontal error 

bars show the experimental results: the left end of each 
error bar corresponds to the case of the exact focusing 
given by the PSF of 2.4 /zm, the right end to the case 
with the PSF of 3.4 /im, and the short vertical bar in 
the middle to the case with the PSF of 2.4 jxm and the 
rotation axis being inclined towards/away from the cam- 
era by 5° [2l[. The two dotted boxes show the range of 
errors in and f2o for the two shots of the experiment. 
Theoretical results are shown by points: the blue squares 
for the PSF of 2.4 /im (exact focusing), the red "+" signs 
for the PSF of 3.4 /im, and the green "x" signs for the 
case with the axial tilt by 5° and the PSF of 2.4 /xm (here 
we referred to as the 5°-tilt case). We observe that the 



the present theory are functions of Rj_ and a (or R z ) 

Writing the experimental data of (y 2 ) and (z 2 ) as 
(y 2 )cx and (z 2 ) cx , and those of the theoretical predic- 
tion as (y 2 )th and (z 2 )th, for each shot of the experiment 
and a are determined by minimizing A/ y 2 (R± , a) = 

((2/ 2 )th - (2/ 2 ) cx ) 2 and AI x ,(R±,a) = ((z 2 ) th - <z 2 ) cx ) 2 
simultaneously. From the obtained value of a, R z (i.e., 
Zq) is readily calculated by Eq. ((7|), and then f2 is de- 
termined by Eq. §5§ using these values of R± , R z , and 
a. Theoretical value of the mode frequency is calculated 
from Eq. ([8]) . We use the same values of the experimental 
parameters as those in Ref. [5|, i.e., ui± = 2ir x 8.3 Hz 
(d± ~ 3.7 jtrai), lo z = 2ir x 5.3 Hz (d z ~ 4.7 /um), and 
a s ~ 5.6 nm for the triplet state of 87 Rb. 

There are two major sources of error in the experimen- 
tal data of the rms radii [2(J. One is an overestimate of 

(y 2 )ex 2 and (z 2 )ex 2 due to defocusing of the imaging cam- 
era. In the experiment, it cannot be completely excluded 
to have the camera defocused by one depth of focus, i.e., 
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FIG. 2: (Color online) Comparison between the rotation ve- 
locity Qo determined in the present analysis and firJ F obtained 
by the TF theory. 



theoretical results for all the three cases are almost on a 
single line irrespective of the way of correcting the data. 

According to Fig. [H the effect of defocusing is large 
enough that our theoretical prediction and experimental 
data agree within the error due to this effect. The effect 
of the inclination of the rotation axial, on the other hand, 
is not so large that it can explain the discrepancy between 
theory and experiment by itself. However, in the fast 
rotation regime, where the frequency shift occurs, this 
effect gives a significant correction in £Iq of ~ 0.002wj_ 
because the cloud is very oblate with the aspect ratio of 
- 0.2. 

It is instructive to compare the rotation velocity de- 
termined here, Oo, with that obtained by the TF theory, 
f2j F , which is invalid when the zero-point energy in the 
z direction is important. From the same data of (y 2 ) e x 
and (z 2 ) cx , we obtain £}q F in the same way as Hq us- 
ing the three-dimensional TF profile whose optical den- 
sity is n opt (y, z) = (4/3)n(0)flj.(l - y 2 / R\ - z 2 /R 2 z f' 2 , 



and a relation between the rotation rate and the as- 
pect ratio given by the TF theory, Oq F /cl>^ = [1 — 
(R z /d z ) 2 /(R±/dx) 2 ] 1/2 , instead of Eq. ©. The abso- 
lute value of £Iq F itself is rather close to that of Oo: dif- 
ference between them is at most < 4 %. However, the 
difference of this amount matters at fast rotation rates, 
where o-> a x shows the rapid change. Thus we compare 
I — flo/u>± and I — SIq/w± in Fig. O Note a striking 
drop at 0.995 < Uq/uj± < 1, which shows the TF theory 
underestimates the rotation velocity in this region. 

Finally we should mention effects of thermal compo- 
nent, which we have neglected in our analysis. Due to 
the low signal to noise ratio in the fast rotation regime, 
it is hard to measure the fraction of the thermal compo- 
nent and we cannot completely exclude its existence [20j . 
Thus it is possible that the thermal component whose 
eigenfrequency is 2oj z drags the oscillation of the conden- 
sate and could act to increase the oscillation frequency 
of the whole cloud towards 2oj z [20] . 

In this work, we have analyzed the axial breathing 
mode frequency using new data of the rms radii obtained 
in a follow-up measurement with attention to remove ef- 
fects of noise [2(|. We have shown that the effects of 
defocus of the imaging camera and an inclination of the 
rotation axis lead to a significant underestimate of es- 
pecially the former one is so large that, within an uncer- 
tainty due to this effect, our theoretical prediction agrees 
with the experimental data. We have also shown that 
TF theory tends to underestimate f2 in the fast rota- 
tion regime, where this theory is no longer valid. Since 
these effects also exist in the experiment of the Tkachenko 
mode Q, they would partly explain a discrepancy be- 
tween theory and experiment for the Tkachenko mode 
frequency in the fast rotation regime [HI, El • 
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